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AN  INVERSE  METHOD  TO  MEASURE  THE  BREATHING  WAVE 
SPEED  IN  A  LIQUID-FILLED  CYLINDRICAL  SHELL 

1.  INTRODUCTION 

The  breathing  (or  bulge)  wave  is  a  fluid-solid  interaction  wave  that  propagates  in  liquid- 
filled  shells  [1],  and  its  wave  speed  is  typically  modeled  as  a  complex,  frequency-dependent 
quantity.  This  wave  can  be  generated  by  applying  axial  excitation  to  one  or  both  ends  of  the 
shell.  Breathing  and  extensional  wave  motion  is  initiated  at  the  location  of  this  excitation,  and 
both  of  these  waves  propagate  longitudinally  in  the  structure.  Additionally,  the  shell  ends 
reflect  such  wave  motion,  and  the  shell  domain  dissipates  some  of  the  energy,  creating  a  spatial 
field  composed  of  partially  standing  and  partially  propagating  waves.  This  wave  motion 
causes  a  change  in  the  cross-sectional  area  of  the  shell,  which  creates  a  pressure  differential  in 
the  fluid.  The  amplitude  of  this  fluid  pressure  is  used  to  measure  the  breathing  wave,  which  is 
strong  at  lower  frequencies  (typically  below  20  Hz  for  most  end- vibrated,  liquid-filled  shells). 
A  previous  technique  to  measure  breathing  wave  speed  [2]  resulted  in  a  single  frequency 
measurement  that  did  not  include  the  loss  (imaginary)  term  of  the  wave.  Additionally,  this 
method  required  sensors  near  the  endcap  of  the  shell,  which  is  an  area  that  is  very  susceptible 
to  acoustic  scattering. 

In  this  report,  an  inverse  method  is  derived  to  measure  the  breathing  wave  speed  of 
liquid-filled  shells.  First,  the  governing  wave  equation  for  particle  motion  is  solved  to  yield  a 
model  of  the  displacement  and  force  in  the  shell.  The  force  and  displacement  equations  are 
then  written  to  correspond  to  two  sensor  transfer  functions.  These  two  equations  and  their 
corresponding  data  are  now  combined,  resulting  in  a  closed-form  value  of  the  extensional  wave 
speed  (including  a  loss  term)  at  each  frequency  where  a  measurement  was  made.  Next,  two 
governing  wave  equations  for  fluid  pressure  are  solved  to  produce  a  model  of  the  pressure  in 
the  shell.  The  pressure  equation  is  then  written  to  correspond  to  the  transfer  functions  between 
five  equally  spaced  hydrophones  (in  the  fluid)  and  a  forward-mounted  accelerometer.  These 
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five  equations  and  their  corresponding  data  are  combined  to  yield  a  closed-form  value  of  the 
breathing  wave  speed  (including  a  loss  term).  The  breathing  wave  speed  measurement  is  to  be 
made  at  low  frequencies  where  the  breathing  wave  is  spatially  coherent  across  the  length  of  the 
shell.  The  amplitude  of  the  breathing  wave  at  higher  frequencies  is  extremely  small  and 
coherent  wavelengths  do  not  exist;  thus  no  high  frequency  measurements  of  this  wave  are 
made.  This  approach  is  intended  for  use  in  the  Axial  Vibration  Test  Facility  (AVTF)  at  the 
Naval  Undersea  Warfare  Center  (NUWC),  Detachment  New  London,  where  liquid-filled 
shells  are  tested,  typically  under  a  tensile  load. 


2.  EXTENSIONAL  WAVE  SPEED  MEASUREMENT 

In  order  to  measure  the  breathing  wave  speed,  it  is  necessary  to  know  the  extensional 
wave  speed.  This  complex,  frequency-dependent  quantity  is  determined  using  the 
measurements  from  the  forward  and  aft  impedance  heads  in  the  AVTF  [3].  Although  these 
measurements  contain  a  breathing  wave  contribution,  it  typically  occurs  only  at  low  frequencies 
and  can  be  discerned  from  the  extensional  wave  effects  with  this  method.  This  behavior  is 
discussed  in  the  experiment  section. 

The  governing  differential  equation  of  the  extensional  wave  is  expressed  on  the  spatial 
domain  as  a  single  wave  equation  with  particle  displacement  as  the  independent  variable: 


dx2 


+  k2U(x,Co)  =  0  , 


(1) 


where  U(x,co)  is  the  temporal  Fourier  transform  of  the  axial  displacement,  x  is  the  spatial 
location  (m),  (0  is  the  frequency  of  excitation  (rad/s),  and  ke  is  the  complex  extensional 
wavenumber  of  the  shell  (rad/m).  It  is  implicit  in  equation  (1)  that  the  fluid-loading  effects  of 
the  extensional  wave  speed  are  incorporated  in  the  axial  model  of  the  shell  with  a  homogeneous 
medium  that  supports  longitudinal  wave  motion.  Although  this  approximation  is  not  sufficient 
to  model  the  radial  motion  of  the  shell,  it  is  an  accurate  model  of  axial  motion  and  its 
corresponding  extensional  wave  propagation.  The  extensional  wavenumber  is  equal  to 
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where  ce  is  the  complex  extensional  (longitudinal)  wave  speed  of  the  shell  (m/s).  The  energy 
attenuation  in  the  shell  is  defined  with  a  structural  damping  law,  and  therefore  the  wave  speed 
is  a  complex  quantity.  The  real  part  of  the  wave  speed  corresponds  to  energy  transmission  and 
the  imaginary  part  corresponds  to  energy  attenuation. 

The  solution  to  equation  (1)  is 

U(x,  co)  =  Ge**x  +  He~ikeX  ,  (3) 

where  G  and  H  are  coefficients  determined  by  the  boundary  conditions  at  the  ends  of  the  shell 
and  i  is  the  square  root  of  -1.  It  is  not  necessary  to  know  these  boundary  conditions  to 
determine  the  longitudinal  wave  speed  of  the  shell.  The  temporal  Fourier  transform  of  the  axial 
force  in  the  shell  is 

F(x,co)  =  ASEX  dU(£~  =  AsExike{GeikeX  -  He~ikeX)  ,  (4) 

where  As  is  the  cross-sectional  area  of  the  shell  (m2)  and  Ex  is  the  effective  longitudinal 
modulus  of  the  shell  (N/m2).  The  known  parameters  in  equations  (3)  and  (4)  are  the  location 
of  the  sensors  ( x )  and  the  frequency  of  excitation  (co).  Although  the  effective  longitudinal 
modulus  (Ex)  is  unknown,  equation  (4)  will  be  rewritten  as  a  ratio  of  forces  permitting  the 
cancellation  of  this  term  and  the  cross-sectional  area  (As).  The  inversion  of  equations  (3)  and 
(4)  at  the  sensor  locations  will  allow  for  a  measurement  of  the  unknown  extensional  wave 
speed  ce.  This  technique  is  described  next. 

The  extensional  wave  measurement  part  of  the  experiment  has  a  forward  and  aft  pair  of 
sensors  to  collect  data  that  are  in  the  form  of  transfer  functions  between  each  pair.  The  position 
of  the  forward  pair  is  defined  as  x  =  0  and  of  the  aft  pair  a sx  =  L,  where  L  is  the  length  of  the 
shell  (m).  The  two  transfer  function  measurements  used  are  the  forward  displacement  divided 
by  the  aft  displacement  and  the  forward  force  divided  by  the  aft  force.  Their  theoretical  form 
can  be  rewritten  using  equations  (3)  and  (4)  as 
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U(0,ai)  _  G±H  _R 
U(L,co )  Ge^1  +  He~'±eL  1 


(5) 


and 


F(0,co)  =  G-H  =  R 

ikpL  zj „ — ikpL 


(6) 


F{L,(0)  GelkeL  -  He~llCel 
where  and  R2  are  transfer  function  data  from  the  experiment.  Equations  (5)  and  (6)  are 
rewritten  as  functions  of  HIG  and  are  set  equal  to  each  other,  yielding 


C0S(k‘L)=% ’ 


(7) 


where  (p  is  a  complex  quantity.  Using  an  angle-sum  relationship  on  the  complex  cosine  term  in 
equation  (7)  and  separating  the  equation  into  real  and  imaginary  parts  results  in 

ReW 


cosh[Im(fcg)L]  = 


cos[Re()L)L] 


(8) 


and 


sinh[Im(fce)L]=  . 


(9) 

sin[Re(fce)L] 

where  Re  denotes  the  real  part  and  Lm  denotes  the  imaginary  part  of  the  corresponding  complex 
quantity. 

Equation  (9)  is  now  squared  and  subtracted  from  the  square  of  equation  (8),  yielding 

{coSh[Im(*e,L]}2  - {stahflm (ML]}2  =  { = 1  '  <10) 

Equation  (10)  can  be  simplified  using  trigonometric  power  relationships  to 
cos[2Re(£e)L]  = 

[Re(0)]2  +  [lm(0)]2  - ^j([ Re(0)]2  +  [lm(0)]2)2  - (2[Re(0)]2  - 2[Im(0)]2  -lj=s  .  (11) 

Note  that  only  a  negative  sign  in  front  of  the  radical  is  used.  The  real  part  of  ke  in  equation 
(1 1)  is  now  solved  for  by 
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Re(U  = 


1  .  ,  .  tin 

— Arccos(^)H - n  even 

2  L  2  L 


— Arc  cos(-5)  +  —  n  odd 

2  L  2  L 


(12) 


where  n  is  a  nonnegative  integer  and  the  capital  A  denotes  the  principal  value  of  the  inverse 
cosine  function.  The  value  of  n  is  determined  from  the  function  s,  which  is  a  cosine  function 
with  respect  to  frequency.  At  zero  frequency,  n  is  0.  Every  time  s  cycles  through  n  radians,  n 
is  increased  by  1.  The  imaginary  part  of  ke  is  determined  by  adding  equations  (8)  and  (9) 
together,  resulting  in 


f  R e(0) 

Im  (0)  1 

[cos[Re(£e)L] 

sin[R e(ke)L] 

(13) 

Now  that  the  real  and  imaginary  parts  of  the  wavenumber  ke  are  known,  the  complex-valued 
extensional  wave  speed  can  be  determined  at  each  frequency  with 


a 


ce  =  Re(ce)  +  ilm(ce)  =  —  . 


(14) 


Note  that  the  extensional  wave  speed  has  been  measured  without  knowing  the  boundary 
conditions  at  x  =  0  and  x  =  L. 


3.  BREATHING  WAVE  SPEED  MEASUREMENT 

The  pressure  field  inside  the  fluid-filled  shell  is  used  to  measure  the  breathing  wave 
speed.  The  dynamic  model  of  this  field  is  derived  from  two  wave  equations  in  the  spatial 
domain,  both  of  which  use  pressure  as  the  independent  variable.  The  first  equation  models  the 
extensional  wave  contribution  and  is  written  as 


d2Pf 


g(x,6))+2p  ,  o 

dx1 


(15) 


where  Pe(x,co)  is  the  temporal  Fourier  transform  of  the  pressure  that  is  generated  by  the 
extensional  wave  and  ke  is  the  extensional  wavenumber  that  was  determined  using  the  method 
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in  section  2.  The  second  wave  equation  models  the  breathing  wave  contribution  and  is  written 
as 


d2Pb(x,co ) 
dx2 


+  kbPb{x,a)  =  0  , 


(16) 


where  PB{x,(d)  is  the  temporal  Fourier  transform  of  the  pressure  that  is  generated  by  the 
breathing  wave  and  kb  is  the  complex  breathing  wavenumber  (rad/m),  which  is  equal  to 


where  c b  is  the  complex  breathing  wave  speed  (m/s).  The  solutions  to  equations  (15)  and  (16) 
are  similar  to  equation  (3),  and  they  can  be  added  together  using  the  principle  of  superposition, 
which  yields  the  total  pressure  in  the  shell  as 

P(x,C0)  =  Pe(x,a)  +  Pb{x,co )  =  AeikeX  +  Be~^x  +  Ce*^  +  De~lkbX  ,  (18) 

where  P(x,a>)  is  the  temporal  Fourier  transform  of  the  pressure  that  is  generated  by  both  the 
extensional  wave  and  the  breathing  wave,  and  A,  B,  C,  and/?  are  wave  propagation 
coefficients  determined  by  the  boundary  conditions.  It  is  now  noted  that  the  pressure  field  at  x 
divided  by  the  forward  accelerometer  is 

Pi*’®)  _  Agikex  +  Be~ikex  +  Ceikbx  +  De~ikbx  (19) 

U 

where  U  is  the  temporal  Fourier  transform  of  the  forward  accelerometer;  A,  B,  C,  and  D  are 
wave  propagation  coefficients;  and  P(x,co)IU  has  units  of  [lPa/g.  It  is  not  necessary  to  know 
the  values  of  A  through  D  to  determine  the  breathing  wave  speed,  as  they  will  be  extracted 
from  the  manipulated  governing  equations.  Five  independent  measurements  of  the  spatial 
pressure  field  are  needed  to  eliminate  these  wave  propagation  coefficients  and  solve  for  the 
breathing  wave  speed.  Without  loss  of  generality,  the  origin  of  the  coordinate  system  is 
defined  as  x  =  0  at  the  middle  (third)  pressure  sensor  (hydrophone).  Equation  (19)  is  written 
to  correspond  to  the  locations  of  the  five  pressure  sensors  as 
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P(  28,(0 )  _  ^-ike28  +  geike28  +  Qg  ikb28  +  pgikb28  _  ^  ^ 
U 

(20) 

P(-8,CO )  _  Ae-ike8  +  Bgike8  +  ^e~llcbs  +  =  £2  , 

U 

(21) 

p®: ^  =  a+b+c+d-s3  , 

U 

(22) 

P(8,CO)  _  AgtM  +  +  Ce*1*5  +De~lkb5  =  S4  , 

U 

(23) 

and 

P(2S,CO)  _  ^gike28  +  Bg-ike28  +  ^gikb28  +  pg-ikb28  _  p  ^ 

U 

(24) 

where  5  is  the  sensor-to-sensor  spacing  (m)  and  through  5$  correspond  to  the  measured 

transfer  function  data. 

Equations  (21)  and  (23)  are  now  added  together  to  yield 

(A  +  B)cos(ke8)  +  (C+  D)cos(kb8)  =  (1/  2)(S2  +  S4)  , 

(25) 

and  equations  (20)  and  (24)  are  added  together  to  produce 

(A+  B)cos(ke28)  +  (C+  D)cos(kb28)  =  (1/2XS!  +  S5)  . 

(26) 

Equation  (22)  is  now  rewritten  with  the  term  (A  +  5)  on  the  left-hand  side  and  substituted  into 

equations  (25)  and  (26),  yielding 

(1/  2)(S2  +  S4)  -  S3  cos  (ke8) 
co  s(kb8)  -  co  s(ke8) 

(27) 

and 

(1  /  2)(5i  +  S3 )  -  S3  cos(fce  25) 
cos(kb2S)  -  cos(ke25) 

(28) 

respectively.  Equations  (27)  and  (28)  are  now  set  equal  to  each  other, 
angle  trigonometric  relationship  to  the  cos(kb28)  term  then  produces 

Applying  a  double 

X  cos2  (kb8)  +  Y  cos (kb8)  +  Z  =  0  , 

(29a) 

where 


7 


TR  11,093 


X  =  (52+54)-253cos(^5)  ,  (29b) 

Y  =  S3 co s(ke8) - (1/ 2)0?!  +S5)  ,  (29c) 

and 

Z  =  [(l/2)(5j  +  5g)  +  iS3]cos(fc£<5)  —  (1/ 2)(S2  "b  ^4) cos (^25)  —  (l/2)(52  +  54)  .  (29d) 

Equation  (29)  is  a  quadratic  form  with  solution 


COS  (kfyS)  = 


-F±VF  -4XZ 


=  V^  , 


where  Y  is  a  complex  quantity.  Equation  (30)  is  similar  to  equation  (7)  and  the  corresponding 


solution  is 


cos[2Re(fy)<5]  = 

[Re(yr)]2  +  [Im(i/r)]2  -  ^([Re(i/0]2  +  [Im(yO]2)2  -  (2[Re(yr)]2  -  2[Im(yr)]2  - 1)  =  r  .  (31) 

Note  that  only  a  negative  sign  in  front  of  the  radical  is  used  in  equation  (31).  However,  both 

the  negative  and  the  positive  signs  in  front  of  the  radical  in  equation  (30)  are  needed.  The  real 

part  of  kb  in  equation  (31)  is  now  solved  for  by 

1  *  ,  ,  .  mit 

— Arc  cos(r)  H — —  m  even 

28  25 

Re(fy)  =  j  ,  (32> 

1  *  .  .  mit  ,, 

— -  Arccos(-r)  H — —  m  odd 

.2 8  25 

where  m  is  a  nonnegative  integer  and  the  capital  A  denotes  the  principal  value  of  the  inverse 
cosine  function.  The  value  of  m  is  determined  from  the  function  r,  which  is  a  cosine  function 
with  respect  to  frequency.  At  zero  frequency,  m  is  0.  Every  time  r  cycles  through  it  radians, 
m  is  increased  by  1.  The  imaginary  part  of  kb  is  determined  in  a  manner  similar  to  equation 
(13),  resulting  in 


t  ft  \-hno  f-  Re(y) _ Im(V/)  ) 

W  5  Sejcos[Re(^)<5]  sin[Re(^)5]J 


(33) 
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Now  that  the  real  and  imaginary  parts  of  the  wavenumber  kb  are  known,  the  complex-valued 
breathing  wave  speed  can  be  determined  at  each  frequency  with 

cb  =  Re(q,)  +  iha(cb)  =  f-  .  (34) 

kb 

Using  this  method  produces  two  wave  speed  measurements  because  of  the  retention  of  the 
positive  and  negative  signs  in  equation  (30).  One  of  the  wave  speeds  is  the  breathing  wave 
speed,  and  the  other  is  the  extensional  wave  speed,  which  was  previously  arrived  at  using  the 
method  in  section  2.  The  extensional  wave  speed  is  typically  at  least  one  order  of  magnitude 
greater  than  the  breathing  wave  speed. 

4.  EXPERIMENT 

Use  of  this  model  with  the  inverse  method  corresponds  to  the  physical  testing 
configuration  in  the  AVTF,  as  shown  in  figure  1.  The  AVTF  has  been  designed  to  provide  a 
simple  procedure  for  testing  long  structures  under  varying  tensions  and  temperatures.  The 
longitudinal  shaker  at  the  forward  end  of  the  structure  provides  axial  excitation  to  the  structure. 
A  rope  attached  to  the  aft  end  of  the  structure  and  a  winch  allows  the  tension  to  be  adjusted.  A 
mass  is  attached  between  the  shell  and  the  rope  to  increase  the  force  levels  and  decrease  the 
acceleration  levels.  This  mass  also  produces  an  impedance  change  at  the  end  of  the  bar  that  is 
sufficiently  large  to  allow  accurate  modeling  of  the  rope  behavior  by  a  spring  and  damper  rather 
than  by  a  continuous  media  expression.  A  rail  from  which  thin  Kevlar  lines  can  be  hung  to 
provide  lateral  support  to  heavy  or  long  test  specimens  runs  the  entire  length  of  the  facility. 
The  entire  unit  is  surrounded  by  an  air-conditioned  PVC  duct  to  permit  temperature-dependent 
testing.  Impedance  heads  are  attached  to  the  forward  and  aft  ends  of  the  structure  to  collect 
data  during  a  test.  Each  impedance  head  consists  of  a  single  axial  force  transducer  and  an 
accelerometer.  Hydrophones  in  the  fluid-filled  shell  collect  additional  data  during  the  test.  A 
load  cell  that  measures  the  tension  on  the  structure  is  located  between  the  rope  and  the  mass. 
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Equipment  used  in  this  experiment  includes  a  Zonic  Corporation  model  1215-10-T- 
ZSP86  hydraulic  shaker,  two  PCB  Piezotronics  model  348A  accelerometers,  two  PCB 
Piezotronics  model  233 A  force  transducers,  a  Lebow  model  31/3  load  cell,  and  a  Cordem 
Corporation  model  1215-RMO  tension  drum. 

In  the  first  part  of  the  experiment,  which  measures  the  extensional  wave  speed,  a 
longitudinally  reinforced  fluid-filled  urethane  shell  containing  five  equally  spaced  hydrophones 
was  placed  in  tension,  as  shown  in  figure  1.  The  shell  had  a  mean  radius  of  0.015  m  and  a 
thickness  of  0.0028  m;  the  internal  fluid  had  a  density  of  760  kg/m3.  The  axial  tension  on  the 
shell  was  890  N  and  the  stressed  length  was  12.0  m.  The  point  mass  had  a  weight  of  13.6  kg. 
The  data  from  the  force  transducers,  accelerometers,  and  hydrophones  were  acquired  in  the 
time  domain  with  a  Hewlett  Packard  (HP)  3562  dynamic  signal  analyzer.  The  analyzer  then 
Fourier  transformed  the  raw  data  to  the  frequency  domain  to  obtain  the  desired  transfer 
functions.  The  test  was  run  with  a  frequency  range  between  3  and  100  Hz. 

Equations  (5)-(14)  were  applied  to  the  experimental  test  data  taken  with  the  force 
transducers  and  accelerometers,  and  the  resulting  extensional  wave  speed  of  the  structure  was 
found.  Figure  2  is  a  plot  of  the  magnitude  and  phase  of  the  ratio  of  the  forward  force  to  the  aft 
force  at  each  measurement  frequency  and  corresponds  to  in  equation  (5).  Figure  3  is  a  plot 
of  the  forward  accelerometer  divided  by  the  aft  accelerometer  versus  frequency  and 
corresponds  to  R2  in  equation  (6).  Figure  4  is  the  calculated  extensional  wave  speed  versus 
frequency.  The  solid  line  corresponds  to  the  computed  wave  speed  and  the  dashed  line  to  the 
ordinary  least  square  (OLS)  straight-line  fit.  The  OLS  fit  was  applied  to  the  data  between  35 
and  100  Hz  to  minimize  the  effect  of  the  breathing  wave  interaction  seen  at  lower  frequencies. 
The  resulting  OLS  fit  was  ce  =  624.7  +  0.8/  (m/s)  for  the  real  part  and  ce  =  69.0  +  0.4/  (m/s) 
for  the  imaginary  part  (where /is  frequency  in  Hertz).  These  same  values  are  used  for  the 
extensional  wave  speed  in  the  following  calculations  in  order  to  determine  the  breathing  wave 
speed. 
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Figure  4.  Extensional  Wave  Speed  Versus  Frequency 
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In  the  second  part  of  the  experiment,  transfer  function  data  of  the  hydrophones  divided 
by  the  forward  accelerometers  were  collected.  The  hydrophones  were  spaced  at  intervals  of 
1.83  m  with  a  distance  of  2.96  m  from  the  forward  end  of  the  shell  to  the  first  hydrophone. 
Figures  5  through  9  are  plots  of  the  hydrophone  pressure  divided  by  the  forward  acceleration 
versus  frequency  and  correspond  to  through  in  equations  (20)-(24),  respectively. 
Equations  (25)-(34)  were  applied  to  the  experimental  test  data,  and  the  resulting  breathing  wave 
speed  of  the  structure  was  calculated.  Figure  10  is  a  plot  of  the  breathing  wave  speed  versus 
frequency.  The  plus  symbols  mark  the  data  from  this  experiment  that  were  obtained  using  the 
method  developed  in  section  3.  The  two  solid  box  symbols  are  data  points  determined  with  the 
phase  angle  between  the  forward  accelerometer  and  the  first  hydrophone.  This  measurement  is 
accomplished  by  first  unwrapping  the  phase  angle  (of  to  a  continuous  function  whose 
ordinate  ranges  from  0  to  1000  degrees.  Next,  regions  that  appear  "straight"  are  identified. 
For  these  experimental  data,  the  regions  are  3  to  5  Hz  and  10  to  25  Hz.  An  OLS  line  was  fit  to 
both  regions,  and  a  single  (propagation)  wave  speed  was  determined  using 

c  =  d™^t,  (35) 

where  c  is  the  corresponding  real  wave  speed  (m/s),  d  is  the  distance  between  sensors  (m),  A / 
is  the  change  of  frequency  (Hz),  and  A 6  is  the  change  of  phase  angle  (degrees).  The  phase 
angle  and  corresponding  straight-line  fits  are  shown  in  figure  11.  The  region  from  3  to  5  Hz 
resulted  in  a  (real)  breathing  wave  speed  of  24.0  m/s  and  the  region  from  10  to  25  Hz  resulted 
in  a  breathing  wave  speed  of  54.8  m/s.  The  shortcomings  of  this  "straight "-line  calculation 
(equation  (35))  are  that  it  cannot  incorporate  extensional  wave  interaction,  requires  single 
direction  wave  propagation,  allows  for  measurement  of  only  the  real  part  of  the  breathing  wave 
speed,  and  is  not  a  very  objective  test  procedure.  Additionally,  it  typically  fails  when  the 
sensor-to-sensor  separation  distance  is  even  moderately  larger  than  that  used  in  this 
experiment.  However,  it  does  provide  a  verification  for  the  real  part  of  the  measured  breathing 
wave  speed. 
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In  figure  10,  it  is  noted  that  the  breathing  wave  for  this  specific  structure  is  spatially 
coherent  from  3  to  13  Hz.  Above  13  Hz,  the  imaginary  part  of  the  measurement  begins  to 
diverge.  The  spatial  coherence  length  is  related  to  the  loss  factor  of  the  circumferential 
modulus  of  the  shell.  A  smaller  circumferential  loss  modulus  would  result  in  a  longer  spatial 
coherence  length  and  a  coherent  measurement  at  higher  frequencies.  Respacing  the 
hydrophone  distances  or  use  of  a  shorter  shell  might  also  result  in  coherent  measurements  at 
higher  frequencies.  Changing  the  extensional  wave  speed  by  ±20  percent  and  recalculating 
the  breathing  wave  speed  produced  a  change  in  the  breathing  wave  speed  of  less  than  1 
percent.  Thus,  the  measurement  method  is  very  insensitive  to  incorrect  extensional  wave 

speeds. 
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Figure  5.  First  Hydrophone  Pressure  Divided  by  Forward  Acceleration  Versus  Frequency 
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Figure  6.  Second  Hydrophone  Pressure  Divided  by  Forward  Acceleration  Versus  Frequency 
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Figure  7.  Third  Hydrophone  Pressure  Divided  by  Forward  Acceleration  Versus  Frequency 
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Figure  8.  Fourth  Hydrophone  Pressure  Divided  by  Forward  Acceleration  Versus  Free 
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Figure  9.  Fifth  Hydrophone  Pressure  Divided  by  Forward  Acceleration  Versus  Freque 
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5.  CONCLUSIONS 

The  breathing  wave  speed  of  a  liquid-filled  shell  can  be  determined  with  five  equally 
spaced  hydrophones.  This  technique  separates  the  spatial  pressure  field  into  two  waves:  an 
extensional  wave  and  a  breathing  wave.  First,  the  extensional  wave  speed  is  determined  using 
force  transducers  and  accelerometers  located  at  the  ends  of  the  shell.  Next,  the  breathing  wave 
speed  is  calculated  using  measurements  from  the  hydrophones  and  the  extensional  wave  speed 
measurement.  The  method  yields  a  complex  value  at  every  frequency  for  which  data  are 
collected  and  for  which  the  breathing  wave  is  spatially  coherent.  This  test  can  be  run  very 
rapidly,  and  the  corresponding  computations  are  not  intensive. 
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